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The magnetoconductance of p-n junctions formed in clean single wall carbon nanotubes is studied 
in the noninteracting electron approximation and perturbatively in electron-electron interaction, in 
the geometry where a magnetic field is along the tube axis. For long junctions the low temper- 
ature magnetoconductance is anomalously large: the relative change in the conductance becomes 
of order unity even when the flux through the tube is much smaller than the flux quantum. The 
magnetoconductance is negative for metallic tubes. For semiconducting and small gap tubes the 
magnetoconductance is nonmonotonic; positive at small and negative at large fields. 

PACS numbers: 75.47.Jn, 73.23.Ad, 73.63.Fg 



Magnetoconductance arises from the orbital and Zee- 
man coupling of electrons to the external magnetic field, 
H. As long as the flux through the crystalline unit 
cell is much smaller than the flux quantum 3>o = hc/e 
magnetotransport may be described in the semiclassical 
approximation pQ ignoring the band structure changes 
due to the presence of H. In most crystals this con- 
dition holds at all experimentally realizable fields. Re- 
cently much attention was focused on magnetotrans- 
port properties of carbon nanotube devices. Because 
of their large radius the magnetic field affects the one- 
dimensional electron spectrum [2 even at relatively weak 
fields. This leads to interesting magnetotransport phe- 
nomena [si m m m [3 

In this paper we show that the magnetic field depen- 
dence of the one-dimensional band structure results in a 
peculiar mechanism of magnetoconductance of p-n junc- 
tions in carbon nanotubes. This mechanism is relevant to 
the magnetoresistance of nominally undoped metallic and 
small gap nanotubes placed on an insulating substrate. 
In this case the long range disorder potential caused by 
charged impurities in the substrate creates p- and n- re- 
gions in the tube, and backscattering of electrons arises 
mainly from the gaps between p- and n- regions, where 
the semiclassical description of electron transport fails. 

We study the magnetoconductance of a p-n junction 
formed in a clean single wall carbon nanotube for mag- 
netic fields parallel to the tube axis. The device is de- 
picted in Fig. [I] a). The p- and n- regions can be formed 
by appropriately biasing the top gates. Such devices were 
recently used to measure [11] thermodynamic properties 
of electron liquid in carbon nanotubes. It is shown be- 
low that for realistic device parameters similar to those 
of Ref. [U the low temperature magnetoconductance be- 
comes of order unity while the flux <1> through the tube 
cross-section is much smaller than 

Before proceeding to detailed calculations let us qual- 
itatively discuss the origin of strong magnetoconduc- 
tance. We choose the x- and y- axes to be respec- 
tively along the tube axis and along its circumference. 
The one-dimensional electron sub-band spectrum is de- 
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FIG. 1: a) Device sketch: a nanotube rests on an insulating 
substrate. The n- and p- regions are created by biasing the 
top gates, b) Band diagram of the device. Tilted solid lines 
represent the bottom of the conduction and the top of the 
valence bands. The width of the center region is 2L. The 
width of the classically forbidden region is 2x Q = 2A/eE. 
Inset: one dimensional spectrum is obtained by dissecting 
the graphene Brillouin zone by the p y = line (dashed line). 
The low energy states lie in the vicinity of K and K' points. 



termined by intersections of the graphene Brillouin zone 
with the k y — ($/<I>o + m)/R lines, see the inset in 
Fig. [I] b). Here R is the tube radius and m an inte- 
ger. The electron spectrum near the K and K' points 
becomes ±a/A 2 + (hvp x ) 2 , where v « 8 x 10 5 m/s is 
the electron velocity in graphene, and A is half the en- 
ergy gap between the valence and conduction sub-bands. 
In the center region between the p- and n- banks the 
external potential is assumed to be approximately lin- 
ear, U(x) = eEx, where e is the electron charge and E 
the electric field. The tilting of electron energy bands in 
the external potential produces a spatial region of width 
2xq = 2A/eE, where electron motion is classically for- 
bidden, see Fig. [!]&). The device conductance is governed 
by the Landau-Zener tunneling across this region. With 
exponential accuracy the tunneling probability may be 
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found in the WKB approximation. The electron momen- 
tum p x along the tube depends on the position as p x = 
^(eEx-e) 2 - A 2 /v. This gives the transmission prob- 
ability T = exp (|lm J p x dx) = exp[— ttA 2 /(hveE)]. 
The rigorous calculation given below shows that the pre- 
exponential factor is equal to unity. The magnetocon- 
ductance arises from the flux dependence of the band 
gap [2], A = A ± hv§/R§o, where A is half-the band 
gap at zero flux and the ± sign corresponds to the differ- 
ent valleys, K and K' . Accounting for electron spin and 
the two valleys and neglecting the Zeeman splitting one 
obtains for the device conductance 
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If the electric field is not too strong, eER <C hv/R, 
the magnetoconductance becomes of order unity while 
the flux is still small, $ <C ^o- I n the case of metal- 
lic tubes, A = 0, the magnetoconductance is negative. 
For semiconducting and small gap tubes the magneto- 
conductance is nonmonotonic; positive at small fields and 
negative at large ones. The conductance maximum is at- 
tained at &rnax ~ ^qAqR/Hv. For semiconducting tubes, 
A ~ hv/3R, this gives <& m ax ~ 3>o/3. For small gap 
tubes the zero flux gap arises only due to curvature ef- 
fects and is rather small, Aq <C hv/R. In this case the 
conductance maximum is achieved at & max <C $o- A 
nonmonotonic magnetoresistance was recently observed 
in Ref. |8] in nanotube devices with a different geometry. 

In the noninter acting electron approximation the mag- 
netoconductance only weakly depends on the tempera- 
ture T as long as the latter is smaller than the Fermi 
energy in the banks. In this regime the energies of 
electrons participating in transport lie in the narrow 
band of width T around the chemical potential. In this 
energy range deviations of electric potential from the 
linear form U(x) ~ eEx are negligible. This results 
in energy-independent transmission coefficient and thus 
temperature-independent conductance. 

In the presence of electron-electron and electron- 
phonon interactions electrons can be transferred between 
the p- and n- regions at finite temperature by thermal 
activation. At T > |A - A | = Hv$/R$ the rate of 
inelastic processes is practically independent of the mag- 
netic field and magnetoconductance arises mainly from 
the tunneling mechanism discussed above. In this regime 
inelastic transfers shunt the tunneling mechanism and 
suppress magnetoconductance. A crude estimate of the 
characteristic temperature T* , above which the magneto- 
conductance suppression becomes significant, can be ob- 
tained in the tunneling regime by equating the activation 
rate, ~ / °° £ exp[-(A + eEx)/T] = ^- exp(-A/T), 
with £i n being the inelastic mean free path, to the tunnel- 
ing rate, ~ exp ^— f^j-^j • According to this estimate the 



noninteracting electron result, Eq. ([!]) provides a good 
description of the conductance for T < eER. 

At zero temperature the finite reflection amplitude at 
the p-n contact leads to the appearance of Friedel oscil- 
lations in the electron density. The additional scattering 
of electrons from the Friedel oscillations in the presence 
of electron-electron interactions gives a correction [12] 
to the noninteracting result for the device conductance, 
Eq. ([I]). This correction is evaluated below to first order 
in electron-electron interaction. It is given by Eq. (15) 
and is plotted in Fig. [2]b). It remains small even if the 
interaction constant, e /Hv is of order unity. 

The device conductance in the noninteracting electron 
approximation, Eq. ([!]), immediately follows from the re- 
sults of Cheianov and Falko [13] for a graphene p-n junc- 
tion that were obtained using transfer matrices. Below 
we present a consideration in terms of wave functions 
that is more convenient for the treatment of the interac- 
tion correction to magnetoconductance. Electron eigen- 
states of energy e obey the Dirac equation, which by an 
appropriate basis choice can be cast in the form 

[U (x) — e — ihva z d x + Aa y ]ifj — 0, 

where o\ are Pauli matrices. Introducing the dimension- 
less coordinate £ = eEx/VeEHv, energy e = e/VeEhv, 
and momenta q = A/\/vheE and fc(£) = (E/(£) — 
e) I V vheE we rewrite this equation as 



fc(£) — id% —iq 
iq fc(f ) + id% 



(2) 



The dimensionless momentum fc(£) changes from — kf at 
£ — > — oo to +kf at £ — > -f oo, where kf ^> 1 is the di- 
mensionless Fermi momentum in the banks. In the center 
region between the banks the coordinate-dependent mo- 
mentum fc(£) is linear in £, fc(£) w £ — e, and the spinor 
amplitudes u and v satisfy the differential equation, 
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(d 2 + a + z 2 )/ = 0. 

Here a = i — q 2 for u : and a = —i — 
introduced the difference coordinate, z = £ — e . The in- 
dependent solutions of this equation are parabolic cylin- 
der functions [14] that can be expressed in terms of the 
confluent hypergeometric function F(a, 7, z), 
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Two linearly independent solutions of Eq. ([2| are 
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where 
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Equation Q conserves the current along the tube axis, 
Ix — ^o~ z 2p. From the form of the current operator it 
is clear that the top/bottom components of the pseudo- 
spinor ip represent the amplitudes of the right- /left- mov- 
ing waves. Thus the scattering states incident from the 
left, ipLi and right, i/jr, can be found by requiring that the 
bottom/top component of the spinor vanish at £ — > ±oo. 
Using Eqs. Q and ([5| and the large distance asymptotics 
of the confluent hypergeometric function [14], 

F( W ^oc)«^(-,)-« + ^ e V- (6) 

we obtain for the scattering states incident from the right 
and left, 
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Here a is given by 
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The z — > ±00 asymptotics of the right-moving wave 
(the top spinor component) in the scattering states ipL 
and ipR are; 



u L (z) 
u R (z) 
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, (9a) 
.(9b) 



From Eq. (9a) one finds the transmission amplitude, 

|2 



to 



lim 



1 



exp 



~~2~ 



This gives the transmission coefficient is % = exp(— 7rg 2 ), 
in agreement with Ref. [13J Expressing q in terms of 
the system parameters using expressions presented in the 
text above Eq.Q, and taking into account the magnetic 
field dependence of the energy gap one obtains the device 
conductance, Eq. ([!]). 



0.1 rt$> 




b) 



0.5 



1.0 



1.5 



FIG. 2: a) Dimensionless electron density as a function of £ for 
q = 0.3, q = 0.6, and q = 0.9. b) The ratio of the correction 
to the transmission coefficient to the interaction constant A, 
Eq. (15), as a function of q. 



Next we evaluate the first interaction correction to 
Eq. g at zero temperature. To first order in interac- 
tion the correction to the device conductance can be ob- 
tained by considering the change in the transmission am- 
plitude for a particle at the Fermi level that arises from 
the additional scattering from the Hartree-Fock poten- 
tial induced by the electron density [12]. The induced 
Hartree-Fock potential has two qualitatively different ef- 
fects on the transmission amplitude: i) By enhancing the 
effective electric field inside the classically forbidden re- 
gion it increases the tunneling amplitude, and ii) It causes 
additional backscattering from the Friedel oscillations in 
the classically allowed region. The analysis below shows 
that repulsive interaction increases the transmission am- 
plitude. 

Using Eqs. ([5| and (|6|) it is easy to show that the spinor 
wave functions in Eq.T4| are normalized to a ^-function 
of the dimensionless energy, 
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Therefore the electron density, upon subtraction of the 
uniform ion background, is 
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In this equation the energy is measured from the Fermi 
level and the electron hole symmetry of the problem was 
used. Plots of electron density, Eq. ( [TT] ), for different 
values of q are presented in Fig. [2] a). The Friedel oscil- 
lations in charge density appear only at finite refection 
amplitude and fall off as l/£ 2 . The extra power of l/£ 
in comparison with the usual one-dimensional Fermi gas 
arises from the linearly growing external potential. Be- 
cause of the fast decay of the oscillations the correction 
to the transmission amplitude is free from infrared diver- 
gences and arises from distances x ~ yJhv/eE. 

Below we consider the case of a short range interaction. 
This should be a reasonable approximation because the 
Fourier transform of the Coulomb interaction depends on 
the transferred momentum only logarithmically. Since 
the characteristic scattering momentum is ^/heE/v the 
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dimensionless interaction constant can be estimated as 
A ~ f- In eE i R , where n is the dielectric constant of 
the substrate. We restrict the consideration to metal- 
lic tubes, for which electron spectra in the presence of 
the flux remain degenerate in the two valleys. Then the 
Hartree-Fock potential is V(£) = 3An(£), where the fac- 
tor 3 = 4 — 1 arises from the spin and valley degeneracy. 

In the presence of the perturbation potential the wave 
incident from the left at the Fermi level, e = 0, can be 
written as i/jl(0 + x(0- The correction, x(£), to the 
wave function satisfies the equation [Hq + V (£)]%(£) — 
—V(£)iI)l(Q w ^h H = £ — id^a z + qa y being the unper- 
turbed Hamiltonian. To first order in perturbation the 
solution of this equation is 



x(0= / dZ'G R (te'Wt')MZ') 



(12) 



where G R (£,£') = (Ho + irj) 1 is the Green's function 
that can be expressed in terms of the spinors in Eq. Q, 
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In order to find the correction to the transmission am- 
plitude one needs only the large distance asymptotics of 
x(0- Therefore only the on-shell part of the Green's 
function will contribute to Eq. (12) at £ — > oo, 
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Since the wave functions ^2 are related to the scattering 
states ipL,R by the unitary transformation Eq. Q the 
sum over i in Eq. (14) can be understood to run over 
i = L, R. Next we notice that for i = L the integral 



in Eq. (14) is purely real. Therefore this term will only 



change the phase of the transmission amplitude, but not 
its modulus. Thus to compute the first correction to the 
transmission amplitude we may replace x(0 by 

The transmission amplitude can be found from the equa- 
tion t = t +lim^ +oo ( ^('_^f • Using Eq. (|9b| we obtain 



for the correction ST to the transmission coefficient 
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It may be evaluated using Eqs. 0, ([5| and (11), and 
is plotted as a function of q in Fig. [2]b). At weak reflec- 
tion, q < 1, the correction to transmission coefficient is 
small even if the interaction constant is of order unity. 
This is because the Friedel oscillations amplitude is pro- 
portional to the reflection amplitude. Thus the noninter- 
acting result, Eq. 0, provides a good description for the 



low temperature conductance of metallic devices. In the 
tunneling regime, q ^> 1, the relative correction to the 
transmission amplitude is expected to be strong because 
of the exponential dependence of the tunneling ampli- 
tude on the effective electric field. Therefore the method 
described above becomes inapplicable. 

In summary, the low temperature magnetoconduc- 
tance of p-n junctions in clean single wall carbon nan- 
otubes was studied in the geometry where the magnetic 
field is along the tube axis. For weak band-tilting field 
E <C hv/eR 2 the magnetoconductance of long, L ^> R, 
junctions becomes of order unity while the flux through 
the tube is much smaller than the flux quantum. In the 
noninteracting electron approximation the device con- 
ductance is given by Eq. 0. The magnetoconductance 
is positive for metallic tubes and nonmonotonic for semi- 
conducting and small gap tubes. The interaction correc- 
tion to the zero temperature magnetoconductance was 
studied to first order in perturbation theory. It arises 
due to the change in the effective electric field in the 
gap between the p- and n- regions and due to the scat- 
tering from the Friedel oscillations. In contrast to the 
one-dimensional Fermi gas, the Friedel oscillations in the 
present geometry fall off as 1/x 2 , which leads to the ab- 
sence of infrared divergence in the correction to the tun- 
neling amplitude. The net correction to the tunneling 
probability is positive. It is given by Eq. ( [T5| ) and is 
plotted in Fig. [2]b). If the reflection coefficient is not 
too strong the noninteracting electron result, Eq. 0, is 
rather accurate even if the coupling constant is of order 
unity. At finite temperatures, in addition to tunneling 
across the classically forbidden region electrons can be 
transferred between p- and n- regions by being promoted 
across the band gap due to inelastic electron-electron and 
electron phonon scattering. The estimate of activation 
transfer rate shows that tunneling processes described 
by Eq. dominate the transport for T < eER. 
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